(19) 



J 




(12) 



EuropdlsclVPPatentamt 

European Patent Office 
Office europ^endes brevets (11) EP 0 978 444 A2 

EUROPAISCHE PATENTANMELDUNG 



(43) Verdffentlichungstag: 

09.02.2000 Patentblatt 2000/06 

(21) Anmeldenummer: 99113611.0 

(22) Anmeldetag: 12.07.1999 



(51) lnt.Cl7: B62D 65/00. B62D 41/00, 
B62D 21/15 



(84) Benannte Vertragsstaaten: 


(71) Anmelder: 


AT BE CH CY DE DK ES R FR GB GR IE IT LI LU 


Vblkswagen Aktiengesellschaft 


MCNLPT SE 


38436 Wolfsburg(DE) 


Benannte Erstreckungsstaaten: 


(72) Erflnder: 


ALLTLVMKROSI 




Scholz, Sven-Peter, Dr. 


(30) Prioritat: 15.04.1999 DE 19917045 


12169 Berlin (DE) 


04.08.1998 DE 19835196 





(54) Verfahren zur Ermittlung der Verformungen und Spannungen einer aus Teilstrukturen 
bestehenden Gesamtstruktur 

(57) Verfahren zur Ermittlung der Verfbrmung einer 
aus Teilstrukturen bestehenden Gesamtstruktur unter 
EinfluB au3erer KrSfte, wobei die Teilstrukturen zumin- 
dest teilweise aus in einem UmformprozeB plastisch 
verformtem Material bestehen und in einem ersten Ver- 
fahrensschritt die Ver^erung von lokalen IVIaterial- 
kennwerten und der lokalen Blechdicken der 
Teilstrukturen in dem jeweiligen UmformprozeB ermitteit 
werden, und wobei in einem weiteren Verfahrensschritt 
die Verformung der Gesamtstruktur unter Einf luB duBe- 
rer Krdfte ermrttelt wird. wobei vor der Ermittlung der 
Verformung der Gesamtstruktur die verdnderten lokalen 
Materialkennwerte sowie die drtlichen Blechdicken der 
Teilstrukturen auf die Gesamtstruktur abgebiklet wer- 
den. 
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Beschreibung 



[0001] Die voriiegende Eiiindung betrrfft ein Verfahren zur Ermittlung der Verformungen und Spannungen einer aus 
Teilstrukturen bestehenden Gesamtstruktur unter EInf Iu3 duGerer Krdfte nach dem Obertegriff des Patentanspruchs 1 . 
5 [0002] In vieien Bereichen der industriellen Technik warden Entwickiungen zundchst einer Computersimulation, bei- 
spielsweise bezuglich der statischen Oder dynamischen Eigenschaften, unterzogen. Ein am Computer entwickeltes 
Modell kann so bezuglich der zu erwartenden physikalischen Eigenschaften beurteilt werden, ohne da6 ein reales 
Modell gebaut und im Strukturversuch getestet werden muBte. 

[0003] Im Bereich der Automobilindustrie ist es Qblich, die Herstellung einzelner Komponenten, insbesondere bel 
10 Blechumfdrmteilen. zunachst zu simuiieren. Meist wird unter Venn/endung eines FEM-Modelles des Ausgangsmaterials 
der UmformprozeB simuliert Dabei kfinnen Verdnderungen der Blechdicke sowie lokale Verfestigungen bei Annahme 
eines geeigneten Materialmodells ermittelt werden. In einem von der Umlbrmsimulation unabhdngigen Simulations- 
schritt wird die Gesamtstruktur des Automobils oder eine aus mehreren Einzelteilen bestehende Teilstruktur des Auto- 
mobils einer Crashsimulation oder Festigkeitsberechnung unterzogen. In diese Simulation flieBen die Ergebnisse der 
75 Umformsimulation nicht ein, es findet allenfalls eine Ubernahme der Geometriedaten statt oder es werden fur belde 
Simulationen die gleichen Geometriedaten der Konstruktion genutzt. 

[0004] Der vorliegenden Erfindung liegt das Problem zugrunde. ein Verfahren der eingangs genanrrten Art zu schaf- 
fen, bei dem die Material- und Geometriedaten des ndchsten Simulationsschritts z.B. des Crash-Simulationsmodells 
ndher an denen der realen Struktur sind. 

20 [0005] Dieses Problem wird durch ein Verfahren nach Patentanspnjch 1 geldst. Erfindungsgemdf3 ist vorgesehen. 
da3 vor der Ermittiung der Verformungen und Spannungen der Gesamtstruktur die verdnderten lokalen Verzerrungs- 
gr03en, aus der insbesondere die Blechdicke abgeleitet wird, sowie die lokalen inneren Zustandsvariablen der Teil- 
strukturen auf die Gesamtstruktur abgebildet werden. Diese Vorgehensweise bietet den Vorteil, da3 der mechanische 
Zustand der umgefbrmten Einzelteile nach der Fertigung als AusgangsgrOBe fur die Berechnung des gesamten 

25 Modells zur Verfugung steht und ausgewertet werden kann. Neben der Anderung der Geometric, insbesondere der 
Blechdicke, durch die Verformung kann ein durch die Fertigungsprozesse verSndertes Werkstoffverhalten berucksich- 
tigt werden. Lokale Verdnderungen in der Umformzonen von Bauteilen mit groBen Ortlichen plastischen Verformungen 
kOnnen die Werkstoffeigenschaften. wie zum Beispiel die FlieBgrenze und den Bauschinger Effekt, wesentlich beein- 
flussen. Diese gednderten Werkstoffeigenschaften werden im allgemeinen durch innere Zustandsvariablen (z.B. ska- 

30 lare und/oder tensorielle VerfestigungsgrOBen) beschrieben und haben EinfluB auf das Verhalten der Gesamtstruktur 
bei einem Belastungsversuch. Verflnderungen der Festigkeit und der plastischen Verfbrmbarkeit einzelner Bauteile wir- 
ken sich unmittelbar auf das Verfbrmungsverhalten der Gesamtstruktur sowie die notwendige Energie zur Erzielung 
einer unzuldssigen Verformung der Gesamtstruktur aus. Bei BerQcksichtigung der Verdnderung der Werkstoffeigen- 
schaften der Einzelteile in den Fertigungsumfbrmprozessen wird eine hdhere Vorhersagegenauigkeit der Simulation 

35 erreicht. Besonders vorteilhaft ist es, wenn die lokalen Materialkennwerte der Teilstrukturen die FlleBgrenze und/oder 
der lokale Verzerrungszustand und/oder weitere lokale innere Zustandsvariable zur Beschreibung des Werkstoffverhal- 
tens und/oder die lokale Verzerrungsgeschwindigkeit sind. Insbesondere die RieBgrenze ist eine entscheidende GrOBe 
zur Bestimmung von Spannungszustdnden, die eine plastische Verformung herbeifuhren. Wird im Materialmodell das 
Verfestigungsverhalten durch eine Spannungs/Dehnungsfunktion beschrieben. so ISBt sich nur durch Angabe des pla- 

40 stischen Verzerrungszustands, insbesondere der plastischen Vergleichsdehnung. und der Verzerrungsgeschwindigkeit 
der Verlauf einer Spannungs/Dehnungsfunktion eines verfestigten Werkstoffes aus der zuvor bekannten Span- 
nungs/Dehnungsfunktion ermitteln. 

[0006] In einer bevorzugten Ausgestattung des Verfahrens ist weiterhin vorgesehen. daB die lokale FlieBgrenze der 
Teilstrukturen anhand des lokalen plastischen Verzerrungszustands. insbesondere der plastischen Vergleichsdehnung, 
45 sowie der lokalen Verzerrungsgeschwindigkeit ermittelt wird. Der plastische Verzerrungszustand ist fOr die Umfbrmung 
von Bedeutung und wird daher in der Umfbmsimulation ermittelt. Die VerknOpfung dieser GrdBe mit einem geeigneten 
Materialmodell liefert so eine einfach zu handhabende skalare GrOBe pro Element. 

[0007] In einer bevorzugten Ausgestaltung des Verfahrens ist weiterhin vorgesehen, daB die Ermittlung der Verande- 
rung der lokalen Materialkennwerte der Teilstrukturen in dem jeweiligen UmformprozeB mit Hilfe eines Finite-Element- 
50 Modells der Teilstrukturen erfolgt, wobei das Finite-Element-Modell aus einzelnen Elementen besteht, die durch Ele- 
mentknoten und Elementfldchenmittelpunkte gekennzeichnet sind. Ein gleichartiges Verfahren kommt bei der Ermitt- 
lung der Verformung der Gesamtstruktur zum Einsatz. Finrte-Element-Verfahren sind gdngige Praxis im Bereich der 
computergestutzten Simulation. 

[0008] Zur Durchfuhrung des Verfahrens kann weiterhin vorgesehen sein. daB zur Abbildung der lokalen Material- 
55 kennwerte sowie der lokalen Blechdicken der Teilstrukturen auf die Gesamtstruktur die Elemente des Finite-Element- 
Modells der Teilstrukturen und des Finite-Element-Modells der Gesamtstruktur einander zugeordnet werden, die den 
kleinsten Abstand voneinander haben, wobei Materialkennwerte des oder der Elemente der Teilstruktur dem zugeord- 
neten Element der Qesanntstruktur zugewiesen werden. Oblichenweise werden in der Umformsimulation Finite-Ele- 
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ment-Modelle mit einer sehr feinen Elementierung benutzt. Fur die Crashsimulation ist es demgegenuber 
wOnschenswert. Bereiche, in denen eine starke Verfbrmung der Gesamtstruktur erwartet wird, feiner zu elementieren 
und Bereiche, die vermutlich kleiner Verformungen unterliegen werden, grOber zu elementieren. 
[0009] Weiterhin kann vorgesehen sein. da3 zur Abbildung der lokalen Materialkennwerte sowie der tokalen Blech- 

5 dicken in der Tellstrukturen auf die Gesamtstruktur die Kanten der Elemente des Rnite-Elementen-Modells der Teil- 
struktur auf das FInite-Element-Modell der Gesamtstruktur projizlert werden. wobel die Materialkennwerte gewichtet 
dem Element der Gesamtstruktur zugeordnet werden. Diese Ausgestaltung des Verfahrens eriaubt eine besonders ein- 
fache Zuordnung der Geometrie des Umformmodells zur Geometrie des Crashsimulationsmodells. 
[001 0] in einer weiteren Ausgestaltung des Verfahrens kann vorgesehen sein. daB die Elemente des Gesamtmodells 

10 entsprechend ihrer plastischen Vergleichsdehnung Massif iziert werden, wobel jeder der so entstandenen Elementklas- 
sen eine Anfangs-FllefBgrenze zugeordnet wird. Die Berechnung der FlieBgrenze anhand eines Materialmodells wird 
durch die Klassifizierung vereinfacht Da das Materialmodell in der Regel aus Versuchen ermittelt wird, bei denen als 
Ergebnis ein Zahlenpaar Spannung/Dehnung anfailt, kflnnen diese Werte recht einfach In eine Klasseneinteilung 
Spannung/Dehnung Qbernommen werden. 

15 [001 1 ] Im fblgenden wird ein Ausfuhrungsbeispiel der vorliegenden Erf indung beschrieben. 

[0012] Fur die Umfbrmungsimulation wird das HlU'sche Materialmodell mit ebener plastischer Orthotropie gewdhlt, 
welches unterschiedliche R-Werte in den Richtungen O"". 45"" und 90"" zur Walzrichtung berOcksichtigt. Die isotrope Ver- 
festigung kann in diesem Model! noch zusdtzlich von der Dehngeschwindigkeit ^ abhdngen. Es werden dazu Span- 
nungs/Dehnungsfunktionen der Art 
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25 benutzt. Da das Verfestigungverhalten dehnratenabhangig sein kann, ist sowohl bei der Umfbnn- als auch bei der 
anschlieBenden Crashsimulation darauf zu achten, daB in der Simulationsrechnung die richtige physikalische Zeit 
berOcksichtigt wird. 

[0013] Die Materlalmodelle nach Gleichung 1 mit einer skalarwertigen inneren Zustandsvariablen erlauben die 
Bestimmung der Verfestigung aus der plastischen Vergleichsdehnung c^p. Diese wird mittels 

30 
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35 

aus dem plastischen Anteil des plastischen Verzerrungsgeschwindigkeitstensors ofj berechnet. Die GrOBe wird fur 
jeden Integrationspunkt uber die Schalendicke ausgewertet. Das Verfahren basiert auf der Annahme einer isotropen 
Verfestigung unter Vernachidssigung der unterschiedlichen Verfestigung uber die Schalendicke durch Waht eines Mit- 
telwerts sowie unter VernachlSssigung des Eigenspannungszustands. Als Modeilparameter verbleiben nur die Scha- 
40 tendicke und die mittiere plastlsche Vergleichsdehnung, die vom FE-Netz der Umfbrmslmulation auf das Modell der 
Crashsimulation ubertragen werden mussen. Die Zuordnung der Modeilparameter auf das Crashnxxiell erfolgt auto- 
matisch mit Hilfe eines Computerprogramms. Dieses bearbeitet folgende Teilaufgaben: 

1. Geometrische Zuordnung der Elemente 
45 2. Berechnung der gewichteten Elementwerte (Blechdickenverteilung) 

3. Klassifizierung der Elemente (plastische Vergleichsdehnung) 

4. Schreiben einer neuen Materialkarte. 

[0014] Als Eingabedaten sind die Ergebnisse aus der Umformsimulation, die zu bearbeitende Eingabedatei der 
50 Crashsimulation sowie eine Werkstoffdatel erforderlich. Aus diesen Daten erstellt das Programm eine ge&iderte Ein- 
gabedatei fur die Crashsimulation. 

[001 5] Die Zuordnung der Elemente des Umformnetzes auf die des Crashnetzes erfolgt unter der Annahme, daB das 
Netz der Umformsimulation feiner ist als das Netz der Crashsimulation. Es wird zundchst von alien Elementmittelpunk- 
ten des Umformnetzes eine Projektion auf das Crashnetz durchgefuhrt. Die Information, welches korrespondierende 
55 Element des Crashnetzes dabei ermittelt wurde. wird fOr jedes Element des Umformnetzes abgespeichert. Anschlie- 
Bend wird das aritiimetische Mittel der Eiementdicken jener Elemente des Umformnetzes gebildet. welche genau 
einem Element des Crahsnetzes zugeordnet werden sollen. Diese gewichtete Elementdicke wird dann in der Eingabe- 
Datel geandert. 
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Patentanspruche 

1. Verfahren zur Ermitllung der Verformung und Spannungen einer aus Teilstrukturen bestehenden GesamtstruWur 
unter Einflu3 duBerer Krdfte. wobei die Teilstrukturen zumindest teilweise aus in einem UmfbrmprozeB plastisch 

5 verformtem Material bestehen und in einem ersten Ver^rensschritt die Verdnderung von lokalen Materialkenn- 
werten und der lokaien Verformung der Teilstrukturen in dem jeweiligen UmformprozeB ermittelt werden, und 
wobei in einem weiteren Verfehrensschritt die Verformung der Gesamtstruktur unter EinfiuB duBerer Krafte ermit- 
telt wird. dadurch gekennzeichnet, daB vor der Ermlttlung der Verformung der Gesamtstruktur die verdnderten 
lokalen Materiaikennwerte sowie die drtlichen Verformungen der Teilstrukturen auf die Gesamtstruktur abgebikiet 

10 werden. 

2. Verfahren nach Anspruch 1, dadurch gekennzelchnet, daB die lokalen Materiaikennwerte der Teilstrukturen die 
FlieBgrenze und/oder der Verzerrungszustand und/oder weitere innere Zustandsvariable zur Beschreibung des 
Werkstoffverhaltens und/oder die Verzerrungsgeschwindigkeit sind. 

15 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB die lokale FlieBgrenze der Teilstrukturen anhand des 
lokalen Verzerrungszustands und/oder weiterer innerer Zustandsvariablen sowie der lokalen Verzerrungsge- 
schwindigkeit ermittelt wird. 

20 4. Verfahren nach einem der AnsprOche 1 bis 3, dadurch gekennzelchnet daB die Ermittlung der Verdnderung der 
lokalen Materiaikennwerte der Teilstrukturen in dem jeweiligen UmformprozeB mit Hilfe eines Finite-Element- 
Modells der Teilstrukturen erfolgt, wobei das Finite-Element-Modell aus einzelnen Elementen besteht. die durch 
Elementknoten und Elementfldchenmittelpunkte gekennzeichnet sind. 

25 5. Verfahren nach einem der AnsprOche 1 bis 4, dadurch gekennzeichnet. daB die Ermittlung der Verformung der 
Gesamtstruktur unter EinfluB duBerer Krafte mit Hilfe eines Finite-Element-Modells der Gesamtstruktur erfolgt. 
wobei das Finite-Element-Modell aus einzelnen Elementen besteht. die durch Elementknoten und Elementf Idchen- 
mittelpunkte gekennzeichnet sind. 

30 6. Verfahren nach einem der AnsprOche 1 bis 5. dadurch gekennzeichnet, daB zur Abblldung der lokalen Materiai- 
kennwerte sowie der lokalen Verformungen der Teilstrukturen auf die Gesamtstruktur die Elemente des Finite-Ele- 
ment-Modells der Teilstrukturen und des Finite-Element-Modells der Gesamtstruktur einander zugeordnet werden, 
die den kleinsten Abstand voneinander haben, wobei Materiaikennwerte des Oder der Elemente der Teilstruktur 
dem zugeordneten Element der Gesamtstruktur zugewiesen werden. 

35 

7. Verfahren nach einem der AnsprOche 1 bis 6. dadurch gekennzeichnet, daB zur Abbildung der lokalen Materiai- 
kennwerte sowie der lokalen Verformung der Teilstrukturen auf die Gesamtstruktur die Kanten der Elemente des 
Finite-Element-Modells der Teilstruktur auf das Finite-Element-Modell der Gesamtstruktur projizlert werden, wobei 
die Materiaikennwerte gewichtet dem Element der Gesamtstruktur zugeordnet werden. 

40 

8. Verfahren nach Anspruch 7, dadurch gekennzeichnet, daB die Elemente des Gesamtmodells entsprechend ihres 
plastischen Verzerrungszustands klassifiziert werden, wobei jeder der so entstandenen Elementklassen eine 
Spannungs/Dehnungsfunktion zugeordnet wird. 
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